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Abstract: This article introduces a novel image 
encryption algorithm based on DNA addition 
combining and coupled two-dimensional piecewise 
nonlinear chaotic map. This algorithm consists of two 
parts. In the first part of the algorithm, a DNA 
sequence matrix is obtained by encoding each color 
component, and is divided into some equal blocks and 
then the generated sequence of Sugeno integral fuzzy 
and the DNA sequence addition operation is used to 
add these blocks. Next, the DNA sequence matrix from 
the previous step is decoded and the complement 
operation to the result of the added matrix is 
performed by using Sugeno fuzzy integral. In the 
second part of the algorithm, the three modified color 
components are encrypted in a coupling fashion in 
such a way to strengthen the cryptosystem security. It 
is observed that the histogram, the correlation and 
avalanche criterion, can satisfy security and 
performance requirements (Avalanche criterion > 
0.49916283). The experimental results obtained for the 
CVG-UGR image databases reveal the fact that the 
proposed algorithm is suitable for practical use to 
protect the security of digital image information over 
the Internet. 

Keywords: Fuzzy integral, Chaotic Maps, Image 
Encrytion.. 

I. INTRODUCTION 

With the rapid development of the networked 
multimedia, communication and propagation 
techniques, transmission of a wide range of digital 
data, from digital images to audio and video files, over 
the internet or through wireless networks has been 
increased. One of the major issues considering the 
digital transmission in this virtual environment is the 
security of digital images, audio, and video files. 
Hence, a growing number of researches have been 
carried out to prevent illegal access to legal data. 
Encryption is usually one good way to ensure high 
security. Due to some intrinsic features of image such 
as bulk data capacity and high correlation among 
pixels, encryption of images is different from that of 
texts. Various digital methods have been proposed for 
image encryption [1-4]. 



The fundamental ideas for image encryption can be 
classified into three major types for spatial domain: the 
value transformation, the position substitution, and 
their combining form. There already exist several 
image encryption methods in spatial domain, among 
which 2D Cellular Automata(CA)-based methods [5, 
6], chaos-based methods[7, 8], the tree structure-based 
methods [9], are most popular. 

Cellular automata [5, 6], is applied into encrypting 
images due to having several strengths. First, large 
number of CA evolution rules has made many 
techniques available for producing a sequence of CA 
data encryption and decrypting images. Second, local 
interactions have made CA suitable for many physical 
streams such as those in image processing and data 
encryption. Third, recursive substitution in CA is 
computationally simple by the sole use of integer 
arithmetic and/or logic operations. An image 
encryption algorithm is recently proposed in [6] using 
recursive cellular automata substitution. The CA 
structure in the proposed algorithms allows an infinite 
number of blocks to be linked together in a chain. 
There are two potential disadvantages associated with 
using this chain. First, the entire chain could be 
garbled by an error anywhere during the transmission, 
and second, the system security could be weakened by 
excessive link dilution. 

Chaos-based cryptographic scheme [10] is an 
efficient encryption first presented in 1989. It has 
many brilliant characteristics different from other 
algorithms such as the sensitive dependence on initial 
conditions, non-periodicity, non-convergence and 
control parameters [11, 12]. The one-dimensional 
chaos system has the advantages of simplicity and high 
security [13, 14], and many studies [7, 15-17] were 
proposed to adopt and improve it. Some of them use 
high-dimensional dynamical systems [14], other 
studies use couple maps [18-20]; however, all research 
is mainly to increase parameters to increase the 
security. As we know, an ideal image encryption 
scheme should be sensitive to the secret key and key 
space should be large enough to make the brute-force 
attack infeasible [21]. 
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Furthermore, it must have some ability to resist 
outer attack from statistical analysis. However, too 
many parameters will influence the implementation 
and will increase computation; for example, the 
proposed scheme in [13] uses eight parameters and 
combines two chaotic maps to shuffle the image pixels 
. The speed and efficiency will also be low using more 
than one encryption scheme to an image. 

In this paper, we have proposed a novel color image 
encryption algorithm based on DNA sequence addition 
combining and coupled two-dimensional piecewise 
nonlinear chaotic map. The architecture of the 
algorithm consists of two stages: The permutation 
stage uses the generated keystream of Sugeno Fuzzy 
Integral (SFI) to change the position and values of 
DNA strings. In the diffusion stage, the three 
components (Red, Green and Blue) of modified image 
are encrypted in a coupling fashion in such a way to 
strengthen the cryptosystem security. 

The rest of this paper is organized as follows. 
Section II describes the fuzzy measure and fuzzy 
integral. In Section III, theory of the DNA sequence 
operation is explained. The proposed coupled two 
dimensional piecewise nonlinear chaotic map is 
discussed in Section IV. In Section V, the proposed 
encryption and decryption algorithm is introduced. 
Simulation results and security analysis are provided in 
Section VI. Finally, the conclusions are drawn in 
Section VII. 

II. FUZZY MEASURE AND FUZZY INTEGRAL 

A. Fuzzy Measure 



Let X = |x 19 x 2 , ,x n | be a finite set associated 

with n attributes on information source space, and 
donate P(X) as the power set X consisting of all 
subsets of X . A set function [22] g: P(x)^[0,l]is 

called a fuzzy measure if the following conditions are 
satisfied: 

1. Boundary conditions: g( ^ = °' g(X) = 1; 

2. Monotonicity: g ^ ~ g< ^ , if A cz B and 
A,BeP(X); 

3. Continuity: lim^ g( A. ) = gClim^ A { ) , if 

{A^jwis an increasing sequence of 
measurable sets. 

Segeno presented the so-called X fuzzy measure 
[22] satisfying the following additional property that 
for all A,B cz X with AnB = § : 

g(A u B) = g(A) + g(B) + Xg(A)g(B) X > -1 (l) 

In general, the value of X can be determined owing 
to the boundary condition of the g x -fuzzy measure. 



This condition reads asg(X) = l, hence, the value of X 
can be found by solving the following: 

ui=n(i+x gi )(2) 

i=l 

Where^e(-l,+oo),^0, and g. =g({x i }) 

is the value of the fuzzy density function. Eq. (2) can 
be calculated by solving the ( n -i)th degree polynomial 
and finding the unique root greater than -1. 

LetA i ={x 1? x 2 , ,x { } . The values of gCA^by 

the fuzzy measure over the corresponding subsets of 
elements can be determined recursively as follows: 

g(A 1 ) = g({x 1 }) = g 1 (3) 

g(Aj ) = gi + g( A M ) + Xg ig (A M ) 1< i < n (4) 

Therefore, in order to obtain X -fuzzy measure, 
there are n parameters g P g 2 >g3>g4 needed to be 
determined in advance. 

B. Sugeno Fuzzy Integral 

Sugeno fuzzy integral of function h computed over 
X with respect to a fuzzy measure g is defined in the 
form: 

E g (h) = i^x[min(h(y i ),g(A i ))] (5) 

Where h(x 1 )<h(x 2 )<...<h(x I1 ), and h(x ) = 0. 

From Eqs. (4) and (5), it is obvious that the calculation 
of the Sugeno fuzzy integral with respect to ?i-fuzzy 
measure requires the knowledge of the fuzzy density g 
and the input value h. 

III. DNA SEQUENCE DESCRIPTION 

In DNA computing, a DNA string is represented by 
a sequence of four basic nucleotides and is usually 
described by letters A( Adenine), T(Thymine), 
G(Guanine), C(Cytosine), where the pairs (A,T) and 
(G,C) are complementary. In the binary, and 1 are 
complementary, so 00 and 11 are complementary, 01 
and 10 are also complementary. In this paper, we use 

C, A, T, G to denote 00, 01, 10, 11, respectively. For 8 
bit grey images, each pixel can be expressed a DNA 
sequence whose length is 4. For instance: If the first 
pixel value of the original image is 220, convert it into 
a binary stream as [11011100], by using the below 
DNA encoding rule to encode the stream, we can get a 
DNA sequence GAGC. Whereas using 00, 01, 10, 11 
to denote C, A, T, G, respectively, to decode the above 
DNA sequence, we can get a binary sequence 
[11011100]. 
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In recent years, algebraic operations such as the coupled two dimensional piecewise nonlinear chaotic 

addition operation, based on the DNA sequence have maps. In order to have good chaotic properties, t and s 

been proposed by researchers [4, 23]. For example, 11 are chosen as 3 and 0.02 respectively [29, 30]. 

+ 10 = 01 or G + T=A. 

V. THE PROPOSED CRYPTOSYSTEM 

IV. COUPLED 2D PIECEWISE NONLINEAR 

CHAOTIC MAP ^' feneration of the initial conditions and parameters 

The proposed image encryption process utilizes a 
A. Nearest-Neighboring Coupled-Map Lattices 256 bit _ long external secret key (K) in p ermu tation and 

A general nearest-neighboring spatiotemporal chaos diffusion stages. This key is divided into 8-bit blocks 
system, also called nearest-neighboring coupled-map (ki) referred to as session keys. The 256-bit external 
lattices (NCML) [24], can be described by: secret-key (K) is given by: 

z n+1 (j) = (l-e)f(z n (j)) + ef(z n (j + l)) (6) K = k p k 2 ,k 3 ,...,k 32 (12) 

Where n = l,2,...., is the time index; j = l,2,...,T is the Based on the analysis presented in section II, there 

T _• 4. 4. - a f • u *.- a ,~ ,. - are four initial inputs of fuzzy integrals and four 

lattice state index; i is a chaotic map, and £ <= (0, 1) is , * r • * 

membership grades for generating pseudo-random 
a coupling constant. The periodic boundary condition num ber bv the SFI 
z n (j + T) = z n (j) is imposed into this system. 



B. Two Dimensional Piecewise Nonlinear Chaotic 
Map 

Two dimensional piecewise nonlinear chaotic map 
with invariant measure are defined as (see [25, 26], for 
the detail): 



For the image of size - WxH, the initial inputs 
h p h 2 ,h 3 ,h 4 are derived as follows, respectively: 



h m =(((k 



((m-l)x8)+l 



v ((m-l)x8)+8 



(D 1 (x n ,a) = x n+1 = 



(7) 



<D 2 (y n ,b 1 ,b 2 ) = 



Yn+1 : 

y„ + i 



4«Xd-x n ) 
l + 4(a 2 -l)x n (l-x n ) 

4b 2 y n (l- yn ) 1 

l + 4(b, 2 -l)y n (l-y n ) 2 '(8) 

: 4h i yAl ~ Vj X. 6 [I,1] (b) 

l + 4(b 2 2 -l)y n (l-y n ) "2 



) + ^(ki))modW)/W 

(13) 

Where the values h i> h 2> h 3> h 4 must be rearranged in 
increasing order [31] and then using these values, this 
scheme generates four membership grades. These 
membership grades can be determined in many 
different ways. Here, we follow the method introduced 
in [32], namely: 



1 



l + h„ 



m=l, 2, 3,4 (14) 



n(x,y) 



IS^^f?^? 1 ?^!?^!^^:^ iS (a similar In diffusion stage, the initial conditions and the 

parameters ( a ,b p b 2 ) of the coupled chaotic map are 

derived as follows: 



calculation has been presented in [27]): 

i Vfr i 



/p 2 



W x (!- X ) (Pi + (l-PiX)) ny/yfX-y) (P 2 + (1-P 2 y)) 

(9) 
According to [28], these maps are ergodic in [0, 1]. 

C. Applying the Coupling Structure to the Two 
Dimensional Piecewise Nonlinear Chaotic Maps 

In the proposed algorithm, we apply the coupling 
structure to the two dimensional piecewise nonlinear 
chaotic maps as follows: 



x (j) = (((k ( , 1)x4+1 e....ek (i . 1)x4+4 )+£(k i ))/w)modi 

i=i 

(15) 

yo(J) = (((k (i+2 ) X 4 + i©--©k (1+2)x4+4 ) + I f (k ; ))/W)modl 

i=l 

(16) 
b m ,,.2 = (((k, 



, m4)x4+25 © ... k (m _ 1)x4+28 ) + X (k; )) / W)modl 



(17) 
a = (y (l) + x (2) + y (2) + x (3) + y (3) + b 1 +b 2 )modl 

(18) 



x n+1 (j) = (1-£)O 1 (x n (j),a j ) + £0> 1 (x n (j + 1),a j ) 

(10) 

y 1Hl (i) = (l-£)<l> 2 (y n (i),bj.bi) + E<l> 2 (y n a+l),b{.bi) D n . 1 1 v . A1 ., 

,+ 1 x B. Design of the Encryption Algorithm 
n =1,2,3, ,L; j = l,2,...,T 



Where £,01,1^ and b 2 are the system parameters; 
x (j) and y (j) are the initial conditions of the 



According to Fig. 1, the proposed encryption 
algorithm can be divided into the following steps: 
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Step 1. First of all, convert the image of size WxH 
into its RGB components, next, transform each 
component into binary matrix, and then carry 
out DNA encoding for the binary matrix 
according to Section III, to obtain encoding 
matrices DNAR, DNAG and DNAB. The size 
of the DNAR, DNAG and DNAB is 
(Wx(Hx4)). 

Step 2. Divide each of matrices DNAR, DNAG and 
DNAB into blocks DNA R (i,j) ,DNA G (i,j), 



DNA B (i,j) 

the size of blocks is 1x4. 



i = 1,2, ,W j = l,2,....,H 



where 



Step 3. Apply external secret key and produce the 
initial inputs h p h 2 ,h 3 ,h 4 and the membership 

grades g 1 ,g 2 ,g 3 ,g 4 according to section V, then, 

generate sequences x n (i), Y n ( j), M(i, 1), N(l, j), 

n = l,2,3,4. 

Step 4. Use the new initial data to iterate SFI once 
using Eqs. (4) and (5): 

Secret Key 

L_ 



Key Generator 



Permutation 



Fuzzy Integral 



Original Image - 



Encoding DNA —► Matrix Blocks —► Add Blocks — ► 



Recombine Blocks & 
Decoding DNA 



Complement j— 



Diffusion 



Coupled Chaotic Map 



Cipher Image 



Fig. 1. Architecture of permutation-diffusion type of proposed image cryptosystem 



M(i,l) = E.modl 



(19) 



Step 10. Sorting X n and Y n in descending order, we 



x m=i,...,4 (i) = ( ARS0nt((E. mod 1) x 10 14 ), (m - 1) x 4)) mod W § et two new sequences X n and Y n . 

(20) Step 11. Let the location value of sequences X , Y 



Remark: ARS(y,x) performs the x-bit right- 
arithmetic-shift operation on the binary sequence y. 

Step 5. Update h p h 2 ,h 3 ,h 4 as follow: 
h m =X m (i)/W, m = 1,2,3,4 (21) 

Remark: The values of n must be rearranged in 
increasing order[31]. 

Step 6. Set i = i + l, If i<W, go to step 4, otherwise, 
Set j=l, and continue the process from step 
Step 7. 

Step 7. Use the new initial data to iterate SFI once 
using Eqs. (4) and (5): 



M(lj) = E.modl 



Y m=1 2 3 4 ( j) = ( ARS(Int ((E . mod 1) x 1 14 ), (m - 1) x 4)) mod W 



(23) 

Step 8. Update h 1? h 2 ,h 3 ,h 4 as follow: 
h m =Y m (j)/W, m = l,2,3,4 (24) 

Step 9. Setj=j+i, If j<=H, go to step 7, otherwise, 
continue the process from step 10. 



be row coordinates and column coordinates of 
DNA R (i,j) , DNA G (i,j), DNA B (i,j), in Other 
words, it can be expressed as DNA ( X » ,y ) , 

DNA G (x' p ,y' q ), DNA B (x p ,y q ). 

Step 12. Add DNA R (i, j) and DNA R (x' p , y' q ) , 

DNA G (i,j)and DNA G (x' p , y' q ) , 

DNA B (i,j)and DNA B (x p ,y q ), according to 

the rules in Section III, obtaining the result as 
blocks BR, BG and BB, respectively. 

Step 13. Carry out the inverse process of the step 1 for 
the decoding matrices BR, BG and BB, we 
will obtain a real value matrices PR, PG and 



PB. 



Step 14. Two sequences M 



(Wxl) 



and N 



(IxH) 



are 



produced by SFI, Performing the multiply 
operation forM, wn andN , we obtain the 

± (Wxl) x ^(lxH) 

matrix Z whose size is WxH .Use the 

following threshold function f ( x ho get a 
binary matrix: 
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[0 ,0<Z(i,j)<0.5 r?4 -x Step 16. Set i=i, F= false, Bitxorl= 0, Bitxor2= 0, 

f (X) = { 1 ,0.5 < Z(i, j) < 1 Bitxor3= 0, C0(1)= 0, C0(2)= and C0(3)= 0. 

Apply external secret key and generate the 
If Z(i,j) = l,P R , P G andP B is complemented, otherwise initial conditions x (i),y (i),x (2),y (2),x (3),y (3) 

it is unchanged. After the complementing operation, and the parameters a ,\ and t> 2 according to 

we get matrices P' R , P' G , and P' B . section IV. 

Step 15. Matrices P' R , P' G , and P' B are transformed into step 17. Step 17: Use the new initial conditions to 
matrices R(wxH)xi>G (WxH)xl andB (WxH)xl , iterate coupled two dimensional piecewise 

respectively. nonlinear chaotic maps once. 



(a) 




(b) 
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Figure 2: Histogram of the original image of Lena in the (a) red (b) green (c) blue, components, Histogram of the encrypted 

image of Lena in the (d) red (e) green (f) blue, components. 



Step 18. In this step, R i ,R i+1 ,G l ,G l+1 ,B i and Bl+1 are \>™ = (x i (2) + y i (2) + b 1 old )mod 1 (33) 

encrypted and the results are store 
matrices q (1), C l+1 (1), Q (2), C l+1 (2), Q ( 
C i+1 (3) using equations (11) and (12). 



encrypted and the results are stored in the b new = (x (3) + v (3) + b old )mod 1 (34) 

matrices C i (l),C i+1 (l),C i (2),C i+1 (2),C i (3) and 2 l n f \Ln r\ ^ 

Bl ^or m=U3 =C m (r)0C m+1 (r) (35) 



Step 19. 
C.(l) = (R.+int(x i (l)xL) + C i . 1 (l) + SBox(Bitxor 1 ))mod256 

(26) 
C. +1 (l) = (R i+1 +int(y.(l) x L) + Q(l) + SBox(Bitxor 1 ))mod 256 

(27) 
C { (2) = (Gj + int^ (2) x L) + C M (2) + SBox(Bitxor 2 ))mod 256 

(28) 
C i+1 (2) = (G i+1 + int^ (2) x L) + Q (2) + SBox(Bitxor 2 ))mod 256 

(29) 
C. (3) = (B { + int^ (3) x L) + C i4 (3) + SBox(Bitxor 3 ))mod 256 

(30) 
C i+1 (3) = (B i+1 + int(yi (3) x L) + Q (3) + SBox(Bitxor 3 ))mod 256 

(31) 

Remark: ^^> ox performs nonlinear transform S-box 
[33] on variable a. 

Step 20. Update a,b 1 ,b 2 ,Bitxor 1 ,Bitxor 2 ,Bitxor 3 and t> 2 as 
follow: 

a new =(x i (l) + y i (l) + a old )mod 1 (32) 



Step 21. Seti = i + 2. Go to step (17) until the elements in 
the R 
exhaust. 



the ^(WxH)xl>^(WxH)xl an( l ^(WxH)xl 



Step 22. If value of the variable F=false, it is set equal 
to true and the operation will start from step 
23, otherwise, the process will be continued 
from step 24. 

Step 23. Set R (W xH)xi>G (WxH)xl ,B (WxH)xl equal to 
reverse of matrices c(i),c(2) and c(3) 
respectively ,next set C0(1)=0, C0(2)=0 and 
C0(3)=0, Bitxorl= 0, Bitxor2= 0, Bitxor3= 0, 
and jump to step 17. 

Step 24. Transform the elements of the matrices 

*^(WxH)xl \*-)i ^CWxH)xl \^)? ^(\VxH)xl \*) 

into C WxH (1), C WxH (2), C WxH (3) and 

output them as color cipher image. 
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It is obvious that the generation of the keystreams 
depends on the keys of cryptosystem, the length of the 
plaintext (L), the width of color image (W) and the 
plaintext through all the color components (R, G, and 
B). Besides, since coupled two-dimensional piecewise 
nonlinear chaotic map has coupling and nonlinear 
structure, the stream cipher of color components 
depend on each other. These features will in turn, 
strengthen the security of cryptosystem. By employing 
this method, we can strongly claim that our proposed 
algorithm does not suffer from the cryptosystem 
weaknesses. 

C. Design of the Decryption Algorithm 

The decryption procedure is similar to that of the 
encryption process except that some steps are followed 
in a reversed order. Therefore, some remarks should be 
considered in the decryption process as follows: 



(b) 




Remarkl: Since the decryption process requires the 
same keystreams as for decryption of the color image, 
the same secret key K = k i> k 2> k 3'-' k 32 should be used 
for decryption. Hence, According to section V, it is 
possible to set the same initial conditions for SFI and 
coupled two dimensional piecewise nonlinear chaotic 
maps. 

Remark2: We can rewrite encryption equations to give 
the pixels' values as follows: 

R.=(C i (l)-int(x i (l)xL)-C i . 1 (l)-SBox(Bitxorl))mod256 

(36) 

G x = (Q (2) - int(x i (2) x L) - C l4 (2) - SBox(Bitxor2))mod 256 

(37) 

B. = (Q (3) - int(Xi (3) x L) - C i4 (3) - SBox(Bitxor3))mod 256 

(38) 



(c) 
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Figure 3: Correlation analysis of two horizontally adjacent pixels: Frames (a), (b) and (c), respectively, show the 
distribution of two horizontally adjacent pixels in the plain image of Lena in the (a) red (b) green (c) blue, components. 
Frames (d), (e) and (f) , respectively, show the distribution of two horizontally adjacent pixels in the encrypted image of 
Lena in the (a) red (b) green (c) blue , components ;obtained using the proposed scheme. 



VI. SECURITY ANALYSIS FOR THE PROPOSED 
ALGORITHM 

A. Histogram 

Image histogram is a very important feature in 
image analysis. From Fig. 2, it is obvious that the 
histograms of the encrypted image are nearly uniform 
and significantly different from the histograms of the 
original image. Hence, it does not provide any clue to 
employ any statistical analysis attack on the encrypted 
image. 

B. Correlation Analysis of Two Adjacent Pixels 

We have analyzed the correlation between two 
vertically adjacent pixels, two horizontally adjacent 
pixels, and two diagonally adjacent pixels in an image. 
5000 pairs of two adjacent (in vertical, horizontal, and 
diagonal direction) pixels from plain-image and 
ciphered image were randomly selected and the 
correlation coefficients were calculated. 



Fig. 3 is the horizontal relevance of adjacent 
elements in image before and after encryption. Fig. 3 
shows significant reduction in relevance of adjacent 
elements. 

C. Avalanche Criterion 

As we know the change of one bit in the plaintext 
should result in theoretically 50% difference in the 
cipher's bits. Hence, for proving the so called 
sensitivity to plaintext, two plain images are generated 
with just one-pixel difference. The bits change rate of 
the cipher obtained by proposed algorithm is 
49.916283%. Hence, the avalanche criterion of is very 
close to the ideal value of 50%. The results of the tests 
comparing the avalanche criterion of the proposed 
algorithm and other encryption algorithms are shown 
in Table I. 
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Table I: Avalanche Test 



Algorithm 



Avalanche 
criterion 



Proposed 
Ref.[19] 
Ref. [20] 



0.49916283 
0.49820101 
0.00000095 



VII. CONCLUSION 

A novel algorithm based on DNA addition 
combining and coupled two-dimensional piecewise 
nonlinear chaotic map has been proposed for the 
permutation-diffusion architecture. The proposed 
algorithm combines good permutation and diffusion 
properties which can be applied to the encryption of 
color images. Combination of image pixels in 
permutation stage based on the generated keystreams 
by Sugeno fuzzy integral is performed. Meanwhile, the 
keystreams used in the diffusion stage is extracted 
from a two-dimensional piecewise nonlinear chaotic 
map. 

The generation of the keystreams depends on the 
keys of cryptosystem, the length of the plaintext (L), 
the width of color image and the plaintext through all 
the color components (R, G, and B). Besides, since 
coupled two-dimensional piecewise nonlinear chaotic 
map has coupling and nonlinear structure, the stream 
cipher of color components depend on each other. 
These features will in turn, strengthen the security of 
cryptosystem. By employing this method, we can 
strongly claim that our proposed algorithm does not 
suffer from the cryptosystem weaknesses. 

Simulation results demonstrate that satisfactory 
performance is achievable in our proposed algorithm. 
Hence, based on the achieved results, we have come to 
a conclusion that the new image encryption algorithm 
with speed and high security can be useful for the real- 
time secure image and video communication 
applications. 
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